A FULLY DISTRIBUTED ALGORITHM FOR THROUGHPUT 
PERFORMANCE IN WIRELESS NETWORKS 

EYJOLFUR I. ASGEIRSSON, MAGNUS M. HALLDORSSON, AND PRADIPTA MITRA 

Abstract. We study link scheduling in wireless networks under stochastic arrival processes of 
packets, and give an algorithm that achieves stability in the physical (SINR) interference model. 
The efhciency of such an algorithm is the fraction of the maximum feasible traffic that the algorithm 
can handle without queues growing indefinitely. Our algorithm achieves two important goals: (i) 
Cn efficiency is independent of the size of the network, and (ii) the algorithm is fully distributed, i.e., 

' ' individual nodes need no information about the overall network topology, not even local information. 
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1. Introduction 



Designing high-performance scheduhng algorithms for wireless networks has become an increas- 

ingly important topic in recent years. Scheduling in a wireless environment is a non-trivial problem, 

I— ( since simultaneous transmissions interfere with each other in complex ways. A two-fold challenge of 

^ appropriately modelling the interference, and then developing algorithms for that model presents 

t/3 itself. Furthermore, in many realistic settings, a centralized controller cannot be assumed, and 

I ^1 algorithms that work in a distributed fashion have to be developed. 

In this work, we are interested in stability and the associated throughput performance of sched- 
'T^ uling algorithms for wireless networks under realistic interference models. We assume that packets 

p^ arrive at potential senders according to a stochastic process, and the goal of an algorithm is to 

\C> schedule these transmissions so that the queues of unscheduled packets at each sender remain 

^^ bounded (in which case, the system is called stable). A rich body of research has been devoted 

; to dealing with this issue in a variety of settings. The seminal work of Tassiulas and Ephremides 

^J^ p3] established that an optimal scheduling policy exists, one that stabilizes the system under all 

pg arrival rates for which stability is potentially possible. In most settings, however, such a "perfect" 

T— I solution is computationally intractable, and additionally a distributed implementation is unlikely. 

^ Hence, the search for efficient and/or distributed algorithms which, if not as good as the optimal 

algorithm, are nevertheless useful. Since these algorithms may not stabilize all feasible arrival 
processes, the concept of efficiency of an algorithm has been introduced, being the fraction of A that 
the algorithm can stabilize, where A is the space of arrival processes that the optimum algorithm 
can stabilize. There have been many approaches to developing such algorithms. A natural step in 
the search for efficient algorithms is to seek maximal solutions. In the context of wireless networks 
this is known as Greedy Maximal Scheduling (GMS) algorithm [10] or Longest Queue First (LQF) 
algorithm [3]. The stability and efficiency of LQF has been investigated extensively [lOl [161 BJ- 
Many other approaches have been proposed as well |181 [TH [22] . 

Most analytic work on wireless networks has been done in graph-based interference models (e.g. 
[18 [llO|[T6l Hll22j^. In these models, wireless links (a link is a sender- receiver pair) that are neighbors 
in a specified link-graph cannot transmit simultaneously. Though interesting in their own right, 
these models are known to over-simplify interference coupling [l7|[20j. As a result, many research 
communities working on wireless networks have increased their focus on the so-called physical model 
or the SINR model. In this model, a transmission is considered successful if the signal received at the 
intended receiver is suitably larger than the cumulative interference due to all other transmissions 
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in the network, plus the ambient noise. In the SINR model, solving the characterization in 
is equivalent to solving the maximum weighted capacity problem, which is known to be NP-hard 
[B] and additionally has no known constant factor approximation algorithm (quite apart from the 
issue of distributed implementation). 

In this paper, we develop an algorithm that is completely distributed, with nodes requiring no 
topological information about the network (not even information about "neighbors"), and achieves 
an efficiency ratio that is independent of the network size (i.e., the number of links). Thus, the 
algorithm is scalable in relation to network size. It can also operate in an asynchronous setting, 
with nodes appearing arbitrarily (as long as the stochastic process meets the required condition). 

We give simulation results which lend credence to the theoretical bounds. Our algorithm is 
extremely robust under fairly high amount of load (achieving efficiency ratios bordering on 0.5). 

The only other work (that we are aware of) on stability of algorithms in the SINR model is the 
recent work by Le et. al. |14j . who analyze the stability properties of the LQF algorithm. The 
authors show that the basic LQF is not efficient, but a variation of it that localizes interference is 
shown to work, with efficiency similar to ours. In a related work, [15] also consider the SINR model, 
but the model there is different (links are always feasible, but have different data rates based on 
the SINR achieved. This sort of problem is rather different from the "combinatorial" situation at 
hand) . 

A distinguishing feature of our distributed implementations is that they require almost no addi- 
tional "infrastructure" . Often distributed algorithms for wireless networks have to assume another 
underlying information infrastructure that can be used to run a localized and/or distributed al- 
gorithm, and that infrastructure, moreover, is not subject to the interference constraints of the 
original network. This is the case with [14], as well as many other works on the topic ([18', [15] 
for example). This is a rather strong assumption, especially in light of the fact that in a wireless 
network, one is usually trying to establish such an infrastructure in the first place. It is interesting 
that we can do without them while obtaining high throughput performance. 

From a technical perspective, we adopt the vocabulary and techniques developed in the context 
of worst-case algorithmic research on the SINR model ([19 l [9l[T2 [ [8j). The concept of "affectance" 
(defined later) developed in some of these works turns out to be quite effective in this context. This 
approach may have further applications in the study of stability of wireless networks. 

The paper is organized as follows. In Section [2| we describe our algorithm and state the main 
result. In Section [3] we present the system model, and discuss related work further in Sec. [4} The 
proof of the stability result is given in Section [5] Finally, in Section [6] we present simulation results. 

2. Algorithm and Result 

The wireless network is modeled as a set L of n links, where each link l^ € L represents a 
potential transmission from a sender s^ to a receiver r„, each a point in a metric space. 

We assume that packets arrive at the sender of each link /„ according to a stochastic process 
with average arrival rate m^. 

The extremely simple and fully distributed algorithm is as follows. 

Our main result is: 

Theorem 1. For all given networks with links on metric spaces, and all sublinear, length-monotone 
power assignments, Reflect achieves an efficiency ratio independent of n. 

3. Some Preliminaries 

The distance between two points x and y is denoted d{x, y). The distance from l^s sender to ly^s 
receiver is denoted d^v = d{su, r^). The length of link /„ is denoted simply by lu = d{su-, ru)- The 
link set is associated with a power assignment P, which is an assignment of a transmission power 
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Algorithm 1 Reflect (Run by each link /„ in the system) 



1: Qu •^ (queue of outstanding packets) 

2: for t ^ 1,2, ... do 

3: Let A be the set of packets that arrive at the beginning of time slot t 

4: Add A to the end of Qu 

5: if Qu is non-empty then 

6: Transmit a packet from Qu with probability 2.5 • rriu 

7: end if 

8: end for 

Pu to be used by the sender of each link lu G L. The signal received at point y from a sender at 
point X with power P^ is Px/d{x, y)" where the constant a > is the path-loss exponent. 

We can now describe the physical or SINR-model of interference. In this model, a receiver Vu 
successfully receives a message from the sender Su if and only if the following condition holds: 

(1) ^"/^" > p , 

S^^,es'\{^„} Pv/d-vu + ^ 

where N is the environmental noise, the constant (3 denotes the minimum SINR (signal-to-interference- 
noise-ratio) required for a message to be successfully received, and S is the set of concurrently 
scheduled links in the same slot (we assume that time is slotted). We say that 5 is SINR-feasible 
(or simply feasible) if M is satisfied for each link in S. 

A power assignment P is length-monotone if P^ > P^ whenever £y > £ui and sublinear if 75- < 7^^ 
whenever i.^ ^ £w This class includes the most interesting and practical power assignments, such 
as uniform power (all links use the same power), linear power (Pu = iu, known to be energy efficient 
in the presence of noise), and mean power {Pu = iu , the assignment that produces maximum 
capacity in this class [81). Let A = -f^^ where £max and £min are, respectively, the maximum and 
minimum lengths of links in L. 

Definition 2. The affectance a^{u) of link /„ caused by another link l^, with a given power as- 
signment P, is the interference of l^ on lu relative to the signal received, or 

a„ (u) = mm <^ 1, c„— • I -— 

where Cu = 13 /(l-PNi^/Pu). 

We need the following assumption. 
Assumption 3. c^ < 2/3 for any link lu- 

This is is fairly reasonable assumption. It simply says that in the absence of other links, the 
transmission succeeds comfortably. The constant 2 is not fundamental; any value greater than 1 
would suffice. 

Since c^ > /? by definition, this implies that 

(2) — < 2 for any two links /„, /„ 

Cu 

The definition of affectance was introduced in [5^ f9] and achieved the form we use in fl3] . When 
clear from the context we drop the superscript P. Also, let a^{v) = 0. Using affectance, Eqn. [ij 
can be rewritten as 

(3) af («) ^ Y. "^(^) ^ 1 ' 
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for all lu & S. 

Signal-strength and robustness. A 6-signal set of links is a set of links where the affectance on 
any link is at most 1/S. A set is SINR-feasible iff it is a 1-signal set. We know: 

Lemma 4 ([7J). Let luAv be links in a q'^-signal set. Then, duv ■ d^u ^ Q^ ■ ^u'^v 

Now we switch to the aspects of this paper related to queueing theory and stochastic processes. 
We first define stability. 

Definition 5. An algorithm stabilizes a network for a particular arrival process if, under that 
arrival process, the average queue size is bounded (at any given time). 

The throughput region is then the set of all possible arrival rate vectors such that there exists 
some scheduling policy that can stabilize the network. 

As proved in [23], the throughput region is characterized by 

A = {A : A ^ (^, for some 4> ^ Co{^)} , 

where ft is the set of all maximal feasible schedules (meaning arrival processes that put weight 
1 on a single maximal feasible set) and Co{Q) is the convex hull of 0,. Note that A and (p are 
n-dimensional vectors and A ^ (/> means each element of A is upper bounded by the corresponding 
one in (p. 

Since fast and/or distributed algorithms might not stabilize all of A, one hopes to achieve a large 
efficiency ratio. 

Definition 6. The efficiency ratio ^ of a scheduling algorithm is 'j = sup{?y : all networks are 
stabilized for all A G ??A}, where rjK = {r/A : A G A}. 

We assume that the arrival process on a link is i.i.d. across time, and different links are indepen- 
dent of each other. 

We will use Mi to denote both maximal feasible sets, and characteristic vectors of said sets (the 
usage being clear from context). For a given efficiency ratio 7, it must hold for all permissible 
arrival rate vectors A that A ^ ^^ rriiMi, where and nrii are weights such that 

(4) ^ mj = 7 . 

i 

It can be easily seen that for any link /«, 

(5) mu= y^ rni<'y . 

4. Related work 

As stated in Thm. [II the efficiency of the algorithm is independent of n (the number of links 
in the system). It is, however, dependent on another network parameter A, the ratio between 
the longest and the shortest link in the system (the proof in the next section contains the exact 
expression). The only comparable work on this model [13] has the same dependence on A (this 
is not explicitly stated in the paper, but can be seen to be necessary). The main discriminating 
feature of our work is that it is distributed in a much stronger sense. The algorithm in [14J can 
be characterized as "localized" , where each link needs to be aware of and have communicated with 
other links in its neighborhood. We have no need for such infrastructure. 

In terms of efficiency ratio, a range of results have been derived in a variety of models. Naturally 
one seeks efficiency of 1 whenever possible [18], but results for efficiency ratio of 1 under certain 
conditions [3], or 2 [3], or g [lOj can be found in the literature. Ratios in terms of certain network 
characteristics are known as well - such as in terms of the degree of the interference graph [2j or 



the local pooling factor [101 114] . In |21| . the abstract SINR model (received signal is a general 
function instead of being length-based) in the context of MIMO networks is studied. An efficiency 
ratio based on a system-specific value ("effective interference number") is derived, with no direct 
comparison with distance-based SINR models. 

For the SINR model, an efficiency ratio that is an "unconditional" constant (independent of both 
n and A) is not known. 

5. Proof of Stability 

We now present a proof of Thm. [T] 

Note that the probability 2.5 • rriu used for link /„ in the algorithm is well-defined, since we claim 
stability with constant efficiency ratio bounded from above by o. 
We first need the following observation. 

Observation 7. For any two links i^ andi^ using a length monotone, sub-linear power assignment, 

P /" 

p fa — 

Proof. If iv 1^ ^u-, it holds by sub-linearity that Pu/^Z ^ Pv/^v- Otherwise, if ^u < ^„, then by 
monotonicity P„ < P^ and by definition of A, ^"/^° < A". D 

The following key lemma shows that no link is affected too much by any single feasible set. 
Lemma 8. Consider a feasible set S and a link ly (not necessarily a member of S). Then, 

(6) Y. "-(^) < '^ • A" , 

hes 

for some constant k. 

Proof. We use the signal strengthening technique of [D]. For this, we decompose the set S to 
[2 • 3"//3]^ sets, each a 3"-signal set. We prove the claim for one such set; since there are only 
constantly many such sets, the overall claim holds (with the appropriate increase in the constant 
factor). Let us reuse the notation S to be such a 3"-signal set. 

Consider the link /„ = (su,r„) G S such that d{ry,ru) is minimum. Also consider the link 
Iw = (swii^w) S 5" such that d{sw,ry) is minimum. Let D = d{ry,ru). We claim that for all links 
Ix = {sx,rx) G S with ix i^ ^w, it holds that 

(7) d[sx.r^)>^B. 

To prove this, assume, for contradiction, that d{sx.,rv) < \D- Then, d{sw,ry) < ^D, by definition 
of Iw. Now, again by the definition of /„, d{rx,ry) > D and d{rw,ry) > D. Thus iw > d{rw,ry) — 
d(rv, Sw) > -J and similarly £x > -f- On the other hand d{sw, Sx) < ^ + ^ = D. Now, dwx ■ dxw < 
{Iw + d{sw, Sx)){^x + d{sw, Sx)) < {iw + D){^x + D) < 9£u>(x, contradicting Lemma[4| 

Now that we have proven Eqn [71 by the triangle inequality, dxu = d{sx,ru) < (i(s^,r„) + 
d{rv,ru) < 3d{ry,Sx) = 3dxv Applying Obs. [TJ we see that 

O-xiv) ^ C„ ^ iji ^ "xu . ZIL <r 9 . qQ P-u-^v ^ o _ Qa\a 

CLx\U) Cu J~^v "'XV ^u ^'o'"u 
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where — < 2 follows from Eqn. l2l Finally, summing over all links in S, 



as{v) = '^ az{v) = aw{v) + ^ az{v) 
iz&S haSMiyj} 

<l + 2-(3A)" Y, «^(^) 

< 1 + 2 • (3A)" , 

where we use aw{v) < 1 by the definition of affectance, and ^^ es\{i \ ^z{u) < 1 since S is feasible 
and lu G -S*. 

This completes the proof setting k = 3"+^. D 

We turn to the proof of Thm. [T| 

Proof. We claim an efficiency ratio of 1/(6kA") where k is the constant from Lemma ^ Thus, it 
is enough to prove stability for all stochastic processes for which the following holds: 



(8) E-^^^ 



6kA" 

Consider the affectance on any link l^ during the execution of the algorithm in a single slot, and 
denote it by a{u). This can be computed as a{u) = ^; Xyav{u), where X^ is a Bernoulli random 
variable which is 1 iff link l^ has a non-empty queue and chooses to transmit during the same slot. 
Now, 

E(a(n)) = ^E(X^,)a„(n) < ^2.5m^,a„(u) 

f'V '•V 



2-5 E 






nii > a^iul 



3orV- Aa 5 2.5kA" 5 
< 2.5 > mjKA" < = — , 

i 

where explanations of the numbered (in)equalities are: 

(1) E(X^) < 2.5mt, by the description of the algorithm. 

(2) By Eqn. Island rearrangement. 

(3) By Lemma p] 

(4) By Eqn. § 

Thus, with probability at least |, a{u) < 1 (by Markov's inequality). Hence, if lu has a non- 
empty queue, with probability at least | • ruu, the queue size decreases. Note that the probability 
can potentially be higher than | • m„, but never smaller. Therefore, this system is at least as 
efficient as the system where in each slot the queue size reduces by 1 with i.i.d. probability exactly 
I ■mu- 

The queue dynamics on a single link become equivalent to the following single server system with 
slotted time and an infinite queue. In this system, at the beginning of each time slot, A packets 
arrive, where A is a random variable on the non-negative integers, with K(A) = m„. At the end of 
each slot, the server processes D packets (or empties the queue), where D is a Bernoulli random 
variable with K{D) = | • rriu > rriu = lE(^). Since the departure process is faster than the arrival 
process, the stability of the queue is guaranteed by basic results in queueing theory [1]. D 

We note that for linear power assignment, the dependence on A in Obs. [7] completely disappears, 
and thus also in Lemma [S] and Thm. [TJ 
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Figure 1. The maximum queue lengths for the distributed algorithms Distr- 
SingleLink and Reflect. The problem instances are based on random topology with 
n = 200, irain = 1, Cax = 20, « = 2.5 and /3 = 1. 



Corollary 9. For all given networks with links on metric spaces, using linear power assignment, 
Reflect achieves an efficiency ratio that is an absolute constant (independent of n and A). 

Remark: In Reflect we have assumed that each link /„ knows m„. In practice, this can be easily 
approximated at time t by min{l, —j^} (where A{t) is the number of packet arrivals up to time t), 
which converges to the right value almost surely. 

5.1. Link partitioning. If link lengths are known beforehand and some pre-processing is allowed, 
the efficiency can be made to have a better dependence on A, specifically, we can achieve an 
efficiency ratio of 1/(6 • 2"Klog A). 

We can partition the link set into a collection of nearly equi- length link sets, i.e., sets where the 
lengths in the set vary by at most a factor of 2. It is easy to show that a link set L can be be 
partitioned into at most log2 A + 1 sets of nearly equi-length links L^ for r = 1 . . . log2 A + 1 where 
Lr contains links of lengths in [2"^'^ • imm, 2*" • ^min)- 

We partition the time slots accordingly, a time slot t is used to schedule links from class Lj. where 
r = {{t- 1) mod (log2 A + 1)) + 1. 

With the partition, the arrival process on a sequence of slots devoted to a single length class 
is equivalent to the setting where all links are nearly equi-length and assuming that Y2i "^i ~ 
1/(6 • 2"«;)rl This partitioning combined with Lemma ^proves the claimed efficiency. 

6. Simulations 

To see how the distributed algorithm Reflect performs, we ran simulations on instances based 
on random topology. The problem instances were created by generating random links in a rectangle 
with side length 100. The length of the links were uniform random variables between ^min and ^max; 
which we set as 1 and 20 respectively. We generated random transmission requests for 100,000 time 
slots while running the algorithms. Our focus was on the behavior of the maximum queue length, 
i.e., the largest number of waiting transmission requests over all the links, which was measured 
every 10,000 rounds. 

Figure [T] shows the results of the distributed algorithm Reflect for random instances with 200 
links after 10,000, 50,000 and 100,000 time slots. The probability scaling factor, p, defines the load 



There are certain technicalities here, since a) Mi n Lt may not be a maximal feasible set in Lt and b) Two sets 



Mi and Mj may have the same "projection" in Lt, i.e. 
can be handled in a straightforward way. 



it is possible that Mi n Lt = Mj n Lt for Mi / Mj . These 



on the system. Thus when p = 1, the system will, in expectation, receive a maximal feasible set in 
every round, while if the scaling factor is zero, no requests are generated. The efficiency ratio of an 
algorithm is then equal to the largest p such that the algorithm is stable. We used a granularity of 
0.01 for values of p between 0.01 and 0.6 and took the average over 10 runs for each value of p. 

The Reflect algorithm used in Figure [T] approximates the arrival rate of requests for each link, as 
mentioned in an earlier remark, instead of assuming that the links know the arrival rate of requests. 

It is interesting to note, in Figure [l| that there is a very sharp threshold where the algorithm 
is no longer stable. As soon as the algorithm becomes unstable, the queue lengths increase very 
rapidly. As expected, the centralized Longest Queue First (LQF) algorithm was more stable than 
the distributed algorithm, managing to keep the maximum queue length below 2 for all values of 
p < 0.6 and not becoming unstable until p > 0.9. We note here again that though our algorithm has 
a lower (though still high) throughput, its main feature is the lack of a requirement for a centralized 
or even localized control. 
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